Crossover from Hydrodynamics to the Kinetic Regime in Confined Nanoflows 
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We present an experimental study of a confined nanoflow, which is generated by a sphere oscillating 
in the proximity of a fiat solid wall in a simple fluid. Varying the oscillation frequency, the confining 
length scale and the fluid mean free path over a broad range provides a detailed map of the flow. We 
use this experimental map to construct a scaling function, which describes the nanoflow in the entire 
parameter space, including both the hydrodynamic and the kinetic regimes. Our scaling function 
unifies previous theories based on the slip boundary condition and the effective viscosity. 



In micron and nanometer scale flows [l|, , the char- 
acteristic dynamic length scale C of the flow approaches 
and is even exceeded by the mean free path of the fluid 
A. This limit is clearly beyond the applicability of the 
Navier-Stokes equations, requiring a rigorous treatment 
using kinetic theory. A less rigorous but widely used ap- 
proach to describe these small scale flows is to extend the 
Newtonian description by imposing a slip boundary con- 
dition on solid walls. This approach is justified as follows. 
Derivation of the Navier-Stokes equations from kinetic 
theory results in the appearance of a Knudsen layer of 
thickness A near the wall Because a fluid element of 
linear dimension ~ A is treated as a mathematical point 
in the hydrodynamic approximation, the velocity at the 
wall becomes u w « A^| z _ Q , with u being the hydrody- 
namic velocity (assumed parallel to the wall) and z being 
the wall normal. Thus, the slip length b, where b ~ A, 
is applied as a convenient empirical parameter to extend 
the Navier-Stokes equations into the kinetic regime. As 
required by macroscopic hydrodynamics, b becomes neg- 
ligible when the Knudsen number, Kn = ^, is small, i.e., 
Kn < 1. 

The above approach comes with some problems. To 
describe some gas flows, for instance, unphysical slip 
lengths, b ^> A, may be required. To alleviate this 
problem, one can assume specular reflections of the gas 
molecules from the wall. Yet, experiments show that 
this assumption is not very accurate for heavier gases 
and untreated surfaces 0, Hf. Worse is the problem when 
the Navier-Stokes solution (with the slip boundary con- 
dition) fails to converge with the prediction of the ki- 
netic theory. A good example to the point is oscillating 
nanoflows [6-8]. Efforts to describe oscillating nanoflows 
using the Navier-Stokes equations in conjunction with a 
slip length agree with experiments only in a range of rel- 
evant parameters A proper kinetic treatment of the 
problem [icl | shows why: the finite relaxation time r of 
the fluid modifies the physics of the flow, resulting in the 
"telegrapher's equation" , which is substantially different 
from the Navier-Stokes equations. 

In this manuscript, we turn our attention to nanometer 
scale confined flows in the limit h < A, where h is the con- 
fining length scale. So far, a group of researchers have 
extended Reynolds' hydrodynamic formulation |ll|, [l2| 



to small scales by imposing the slip boundary condition 
[l3l - [l7} — as described above. Others, coming from ki- 
netic theory, have developed the concept of the effective 
viscosity, which typically depends upon a properly de- 
fined Knudsen number [3 El. There is no question 
that both approaches must agree for the same flow pa- 
rameter space. Here, we present an experimental study of 
nanometer scale confined flows covering a broad range of 
parameters — including gap h, pressure p, and frequency 
— along with a scaling theory. Our scaling function 
describes the physical behavior of the flow in the entire 
parameter space, capturing the transition from hydrody- 
namics to the kinetic regime accurately. 

We study the oscillatory hydrodynamic response of a 
sphere in the proximity of a solid surface. Our experi- 
mental device is a micron-scale silica sphere with radius R 
glued to the end of a microcantilever of linear dimensions 
I x w x t. We have employed both the fundamental and 
first harmonic flcxural modes of a soft cantilever (CI), 
and the fundamental flexural mode of a shorter stiffcr 
cantilever (C2). Figure HJa) depicts optical measure- 
ments and finite element method (FEM) simulations of 
the mechanical modes of CI. For each device and mode, 
we first extract the the intrinsic quality factor Qi and 
resonance frequency ^ in UHV away from any surfaces. 
The modal mass m e is determined from the resonance 
frequency shifts before and after the sphere is attached 
to the cantilever. These parameters are listed in Table Q] 

Once the mechanical mode is characterized, we change 
the flow parameters while optically monitoring the dis- 
sipation and the resonance frequency of the mode. In 
particular, we continuously vary two parameters for each 
mode as follows, i) We change the gap h (shortest dis- 
tance) between the sphere and a flat solid (Silicon) sur- 
face. At small gaps [h < 200 nm), we drive the cantilever 
to achieve 'intermittent contact' between the sphere and 
the solid, and determine the gap from the amplitude. 
For large gaps, h is extracted from a calibrated linear 
motion stage, ii) We vary the surrounding pressure p 
by admitting dry N2 into the chamber. These provide 
a two-dimensional parametric map of the dimension- 
less dissipation and the (angular) resonance frequency: 
Qm~ l = Qm (h,p) and w m = oj m (h,p). Before present- 
ing the data, we show in Fig. QJb) that 1/Q m measured 
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FIG. 1. (a) Measured (symbols) and simulated (inset and solid lines) mechanical modes of the sphere-cantilever device. In 
the first harmonic mode (blue), FEM simulations suggest that the node appears at the position where the sphere is attached 
to the cantilever and that the sphere undergoes small rotational oscillations about an axis (parallel to the y-axis) through the 
node, (b) Measured dimensionless dissipation 1/Q m as a function of gap h, obtained from the driven frequency response, and 
thermal oscillations of the device, (c) Dimensionless fluidic dissipation 1 /Q / as a function of gap h at fixed pressures p. Solid 
line segments show asymptotic values 1/Q/oo from viscous theory. Inset is a semi-logarithmic plot of the same data. The h 
error bars for h < 200 nm are due to roughness and contact uncertainty. Otherwise, the error comes from the linear stage. The 
error bars in 1/Qf are smaller than the symbol sizes, (d) 1/Qf measured as a function of pressure p with the gap fixed. Solid 
line is from viscous theory and the dotted line is from molecular theory. All data in this figure are obtained from device CI. 



by linearly driving the resonator and by monitoring its 
thermal fluctuations agree closely, with a typical discrep- 
ancy less than 1%. The maximum amplitudes in driven 
and thermal measurements remain ~ 1 nm and ~ 0.01 
nm, respectively. By properly subtracting the intrinsic 
dissipation from the measured dissipation, one can ob- 
tain the fluidic dissipation: 1/Q/ = 1/Q m — 1/Qi- Fig- 



ure []Jc) and (d) show the Q/~ 



Qj (h,p) data set 



for the 13.7 kHz mode in double-logarithmic plots against 
gap h and pressure p, respectively In Fig. [IJc), the gap 
is varied in the range 10~ 8 m < h < 10 -4 m with the 
pressure held at p = 100, 300 and 1000 Torr. Conversely, 
in Fig. Hid), the pressure is swept continuously in the 
range 10~ 2 Torr < p < 10 3 Torr, while the gap is fixed at 
h = 0.1, 0.2, 1 and 20 /im. The inset in Fig.[ljc) is a semi- 
logarithmic plot, showing the characteristic saturation of 
Qf~ vs. h (see discussion on 1/Q/oo below). The ac- 
companying mode frequency, oj m = u> m {h,p), show very 
little variation (less than 0.1%) in this parameter space. 

Several important preliminary observations can be 
made from the data of Fig. [TJc) and (d). For a sphere 
oscillating at frequency^- in an unbounded fluid at the 
viscous limit lot <C 1 [6|7 1 1 0I| . the dimensionless dissipa- 
tion can be written as 1201 



1 



Qfc 



6irp,R 
mco 



1 R 

1 + T 



(i) 



where R is the radius and m is the mass of the sphere, 
p is the dynamic viscosity of the fluid, and 6 = W -£j 



TABLE I. Mechanical properties of the measured devices. 



Device 


Mode 


I X W X t 


R i£i. 

(Ami) (kHz) 


Qi 


m e 
(kg) 




CI 


1 


230 x 40 x 3 


35 13.7 


12xlO a 


5x10" 


10 


CI 


2 


230 x 40 x 3 


35 45.8 


3.4xl0 3 


16x10 


-10 


C2 


1 


125 x 35 x 4 


21.5 122.4 


6.8xl0 :i 


1x10" 


10 



is the viscous boundary layer thickness, p = pis, where 
p is the density and v is the kinematic viscosity of the 
fluid. The fluidic dissipation from the rectan gula r can- 



tilever can also be found, albeit numerically [21|. The 
solid line segments in Fig. QJc) and the solid curve in 
Fig. [TJd) show the 1/Q/oo predictions of viscous theory 
at large gaps, h — > oo. In these calculations, the inde- 
pendent contributions to dissipation from the sphere and 
the cantilever arc simply added. The velocity field of 
an oscillating sphere-cantilever system should be differ- 
ent from that obtained by adding the individual velocity 
fields of a sphere and a cantilever. Regardless, the agree- 
ment between experiment and calculations in Fig [TJc) 
and (d) is satisfactor 
of molecular theory 



22] 



In Fig. Q] (d), the prediction 
is also shown. 
When a wall is placed in the proximity of an oscillating 
sphere, the entire velocity field (not just the field in the 
gap) will be modified substantially. Regardless, the dissi- 
pation caused by the squeezing of fluid in the gap can be 
conveniently studied by subtracting the dissipation in an 
infinite fluid, i.e., 1/Q/j = 1/Q/ — 1/Q/oo- Subtracting 
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FIG. 2. Gap dependent dimensionless dissipation l/Qh as a function of gap h at fixed pressure, (a), (b) Fundamental modes 
of devices CI and C2, respectively, (c) First harmonic mode of device CI. Solid lines in (a)-(c) are fits to Eq. with a = 0.5 
and j3 = 1.6, multiplied by a fitting factor of C ~ 0.23 ± 0.11. The deviation from the solid line in (c) is possibly due to the 
additional rotational motion of the sphere. The dashed line in (c) is the improved fit with the added rotational dissipation [24l j. 
The noise in all the data in (a)-(c) increases for h > 10 4 nm due to the subtraction of 1/Q/oo- The representative error bars 
are found by an analysis of the noise at the tails (h > 10 4 ) and become smaller than symbols for h < 10 3 nm. 



the experimental /i-indcpcndcnt l/Q/oo asymptotes in 
Fig.QJc) from the 1/Qf data results in the dimensionless 
gap-dependent dissipation l/Qh in Fig. HJa). Figure [2] 
depicts similarly-obtained 1 / Qh for three different modes 
at multiple pressures as the gap is varied. Solid lines are 
fits to theory (see below) . A first pass analysis of the data 
can be provided based upon the dimensionless Knudsen 
number, Kn/j = ^. When Km, <C 1, l/Qh oc 1/h and 
can be approximated as 23l425j : 
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At the opposite limit of Knj, ^> 1, the dimensionless dis- 
sipation saturates. Between these two limits, there is 
a well-defined transition from the hydrodynamic to the 
kinetic regime. 

We now provide a theoretical background for the ob- 
served transition. Since l/Qh —> as h —> oo, we can 
write a general relation 
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The scaling function f({xi}), which is analytic in the 
limit {xi} — > 0, depends on various dimensionless vari- 
ables pertaining to different dynamic regimes. It is clear 
that the first few {xi} are the familiar Knudsen num- 
bers based on appropriate linear dimensions characteriz- 
ing the system: Ke/, = ^ , Knj = j , Knjj = and so on. 

The last dimensionless variable, = R x /~^ = \/— = 
Re^ , can be regarded as a Reynolds number based on the 
velocity U = uR/2. In the limit Kn^ — > and Kcs — > 0, 
Taylor expansion gives 
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where 



a { ^Re s (5) 
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The relative magnitudes of the a\ and the dimensionless 
parameters {a;,} determine the physics of the flow. By 
varying {xi} over a broad range, one can extract the 
magnitudes of af^ from experiment. 

To gain more insight into the proposed expansion in 
Eq. let us consider the limits. When h — > oo 

(Kn/j — > 0), dimensionless dissipation due to squeezing 
disappears, l/Qh 0. This suggests that the first or- 
der term in the Taylor expansion in Eq. ([5]) should not 
strongly depend on the other Knudsen numbers, Kn^, 
KnR and so on. In the limit of small h (Kn/, S> 1), mo- 
mentum transfer is dominated by the ballistic impact of 
the molecules emitted from the stationary plate incident 
on the moving sphere. The contribution of intermolecu- 
lar collisions can be neglected. If the thermal molecular 
velocity uth is large, the dimension of the gap h must 
disappear from the expression for dissipation. 

Keeping a finite number of terms in Eq. (j4]), one can 
only hope to find an approximation for the scaling func- 
tion f({xi}) valid in the limit xi — > 0. To obtain an 
expression valid in the entire range of {xi} variation, one 
has to keep infinitely many terms. This can be achieved 
by recasting the scaling function in Eq. (J3j) into a ratio 
of low-order polynomials with unknown coefficients to be 
determined experimentally. The resulting expression 
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can be perceived as the simplest Padc approximant, 
which should describe experiments in a broad paramc- 

in) 

ter range. The constants a and j3 are related to a\ . 
It is interesting to note that, in this choice, the term 
of linear order 0(Res) disappears due to the subtrac- 
tion, 1/Qh = l/Qf — 1/Q/oo- However, the higher or- 
der term ©(Kn^Re^) survives. In the small-gap limit 
Knf, = ^ ^> 1, one obtains as prescribed 
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Returning to Fig. [3J we now describe how the fits to 
the experimental data are obtained based upon the above 
scaling form. The device parameters m = m e , uj, and R 
are experimental constants. The fluid parameters are all 
assumed to be independent of h, but may depend on p: 
and fj. is independent of p. The 
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very same constants a and f3 in the scaling function in 
Eq. ([7|) must uniquely fit all data sets — regardless of 
pressure, frequency, mode and so on. Indeed, we can fit 
all our data with a = 0.5 and j3 = 1.6, found by iteration. 
Any small changes in a and j3 cause the curves in Fig. [2] 
to shift along the /i-axis, making the fits unacceptable. 
The fits can be improved along the 1/Qh-axis by mul- 
tiplying with fitting factors of C ph 0.33, 0.20 and 0.16 
for the 13.7 kHz, 122.4 kHz and the 45.8 kHz modes, 
respectively, resulting in the solid curves in Fig. [2] To 
within our experimental accuracy, however, C remains a 
constant as C ps 0.23 ± 0.11 for all our devices, and may 
be needed due to non- idealities in geometry (e.g., the 
cantilever and epoxy above the sphere), inaccuracies in 
determining m e (especially for the first harmonic mode) 
and deviations from normal relative motion (see below). 
Deeper physical factors — such as the non-trivial effects 
of the subtraction of the 1/Qfoo tails and unsteady cor- 
rections to Eq. ((5]) — cannot be ruled out, and may give 
rise to the small deviations in C from device to device. 

The fit in Fig. [3Jc) (solid curve) deviates from the 
data for 10 2 nm < h < 5 x 10 4 nm. FEM simula- 
tions for this mode suggest that the sphere undergoes 
rotational motion — with the displacement of its closest 
point to the wall being in the direction 0.98x + 0.2z [see 
Fig. da)] • Then, the dissipation comes from shearing the 
fluid in the gap as well as from squeezing it. For shear, 
= 4 i5mu x m (f ) as PP ose d to the expression 
in Eq. ^ for squeezing. The dashed line in Fig. [2Jc) 
is the fit found by naively adding these two forms in the 
ratio of the FEM motional amplitudes, and by keeping 
the scaling function exactly the same. Because the effect 
remains small, the 1/h dependence of the dissipation can 
be assumed prevalent for all devices considered here and 
in the literature. 

Having fit individual data traces, we can collapse all 
our data as shown in Fig. [31 The collapse is obtained 
by removing the trivial effects of the device size and frc- 
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FIG. 3. Collapse of all experimental data from this work. 
Here, / is the scaling function defined in Eq. (0 and the 
same Cs are used as above. 



quency from the data as well as the more profound effects 
of the scaling function f({xi}). The plotted dimension- 
less quantity, -q- x 6 ™ U R x A can be regarded as the di- 
mensionless size- and frequency-independent dissipation, 
in which the kinetic effects have been deconvoluted. It 
therefore shows the hydrodynamic R/h dependence at all 
length scales studied here. 

Finally, our results can be interpreted as follows. In 
the hydrodynamic limit (h 3> A), this problem is de- 
scribed by Eq. ([2]), where the viscosity fi is dominated 
by intermolecular collisions, fi ~ pu t h\, with a relax- 
ation timescale ~ \/u t h- To gain insight into the kinetic 
limit (h <C A), one can simply write the shear stress 



on the sphere as cr 



pu t h 



It is easy to see that 



d ~ puthh^jr ~ puthh^, where ^ is the velocity gradi- 
ent. This result can be interpreted as the appearance of 
an effective viscosity, p e ff ~ P u thh, due to an effective 



mean free path, A 



eff 



h. Substituting p e ff m to the 



hydrodynamic solution simply results in -A- 



Qh ' mu ' 

consistent with Eq. ([5]). Thus, in principle, one may 
justify an attempt to reach the kinetic regime by us- 
ing the Navier-Stokes equations, but combined with ef- 
fective (and sometimes frequency-dependent) viscosities, 
slip lengths and so on. 

In this manuscript, we have presented experimental 
data on confined nanoflows covering a broad range of 
flow parameters. Our simple scaling theory describes 
experiments in the entire parameter range — without 
explicitly employing an effective viscosity and/or slip 
length. To conclude, we stress that the dimensionless 
Wcisscnberg numbers here remain small, Wi = lot <C 1. 
Since the appearance of frequency in effective viscosity 



5 



essentially leads to a modification of the equations of mo- 
tion 0, [l(| , generalization of confined nanoflows to the 
interval Wi ^> 1 will require further experimental and 
theoretical work. 
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